






















AD A 1 O 4 6 0 7 


UNCLASSIFIED j rllf ( JfTf /j &) 

SECURITY CLASSIFICATION OF THIS PAGE fPh.n Dmta Ent,r,d) | ■ V P 1 (/ J m- 

DPDftDT noniUCklTATlftlJ DACC VlH V fli 

lH"' READ INSTRUCTIONS’^. — 


BEFORE COMPLETING FORM 

1 REPORT NUMBER 2. GOVT ACCESSION NO. 

r L. M>- 

3. RECIPIENT'S CATALOG NUMBER 

>7 

1 Jj^ftTLE (and SublllleJ 

Some Concepts of Negative Dependence r 

\.. . 

U ' - _ r . - . . .. 

5. TYPE OF REPORT « PERIOD COVERED 

Technical { ; ' 

S. PERFORMING ORG. REPORT NUMBER 

A Au ThORC»> / 

Henry W. Block 

Thomas H. Savits J r 

I Moshe Shaked \ 

e. contract or grant numbero> 

"■ NJTO014-76-C-0839J ^ 

-m 042-357 

9 PERFORMING ORGANIZATION NAME Al D ADDRESS 

University of Pittsburgh 

Department of Statistics 

Pittsburgh, PA 15260 

to. PROGRAM ELEMENT. PROJECT, TASK 
AREA A WORK U4UXJiUU*CRSt 

;. m 


1 I. CONTROLLING OFFICE NAME AND ADDRESS 

Office of Naval Research (411(SP)) 

800 N. Quincy Street 

Arlington, VA 22217_ 



REPORT OATE 

MarepsftjSO ! 

NUMBER OF PAGES 


36 


14. MONITORING AGENCY NAME h AODRESSf// different from Controlling Office) [ IS. 


SECURITY CLASS, (of thte report) 


UNCLASSIFIED 


15«. DECLASSIFICATION/DOWNGRADING 
SCHEME 


16. DISTRIBUTION STATEMENT (of thla Report) 

Approved for public release; distribution unlimited. 


17. DISTRIBUTION STATEMENT (of the abstract entered In Block 20, II different from Report) 


SEP 2 8 198i ■ 


oZ 

18. supplementary notes 

19. KEY WORDS ("Continue on reveree aide If neceeeary and Identify by block number) 

Multinomial Multivariate Hypergeometric Negative and Positive 

Dirichlet Multivariate Normal Dependence 

Inequalities PF Densities New Better Than Used 

Stochastic Ordering Total Positivity 

1 _-I_ 


20. A&^Yr ACT (Continue on reveree aide It neceeeary and Identity by block number) 

"'ihe theory of positive dependence notions cannot yield useful results for 
some widely used distributions such as the multinomial, Dirichlet and the 
multivariate hypergeometric. Some conditions of negative dependence that are 
satisfied by these distributions and which have practical meaning are Intro¬ 
duced. Preservation results for some of these concepts are derived. Useful 
inequalities for some widely used distributions are obtained. Results of 
Mallows (1969) that apply to the multinomial distributions are extended to 
more distributions. Examdes are listed. 


DD | j2J M 7 j 1473 COITION OF 1 NOV 65 IS OBSOLETE A* UNCLASSIFIED 

c /M mm. ip. nu. aaai ■ L ^ i ■ — 


SECURITY CLASSIFICATION OF THIS P AGE Data Bnterad) 


L i O'/ '- r / 











^Department of Mathematics and Statistics, University of Pittsburgh, 
Pittsburgh, Pennsylvania 15260. 

2 

Department of Mathematics, Indiana University, Bloomington, Indiana 47505. 
^Supported by ONR Contract N00014-76-C-0839. 

^Supported by NSF Grant MCS77-01458 and ONR Contract N00014-76-C-0839. 



ft 








Abstract 


The theory of positive dependence notions cannot yield useful 
results for some widely used distributions such as the multinomial, 
Dirichlet and the multivariate hypergeometric. Some conditions of 
negative dependence that are satisfied by these distributions and which 
have practical meaning are introduced. Preservation results for some 
of these concepts are derived. Useful inequalities for some widely 
used distributions are obtained. Results of Mallows (1969) that apply 
to the multinomial distributions are extended to more distributions. 
Examples are listed. 
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1. Introduction . 

Concepts of positive dependence of sets of random variables (rv's) 
have received a lot of attention recently. Their study was found to 
yield a better understanding of the structure of some widely used multi¬ 
variate distribution functions (df's). In addition to this, various 
useful Inequalities were obtained with applications in many areas of 
probability and statistics. Barlow and Proschan (1975), Ch. 5, Include 
a review of most of the work done prior to 1972. A list of more recent 
references can be found in Ahmed et al. (1978). 

On the other hand notions of negative dependence have received very 
little attention in the literature. Some negative dependence analogs of 
positive dependence concepts have been mentioned by some authors (Lehmann 
(1966), Brindley and Thompson (1972), Dykstra, Hewett and Thompson (1973) 
and Shaked (1977) among others). In the bivariate setting the random 
vector (T^^) * 6 usually said to satisfy some negative dependence 

condition if (T^,-^) satisfies the analogous positive dependence 
condition. However, this method of formulation cannot apply to higher 
dimensions. To the best of our knowledge Lehmann (1966) in the bivariate 
setting and Mallows (1968) in the multivariate .setting came the nearest 
to a systematic study of negative dependence concepts; Mallows' discussion, 
however, is restricted to the multinomial distribution. The development 
below is in the spirit of Mallows; his results are special cases of ours. 

While Che first draft of this paper was being written two related 
works were brought to our attention. The first work by Ebrahimi and Ohosh 
(1980) discusses some negative dependence analogs of well known positive 
dependence concepts. Some of our definitions overlap those of Ebrahimi and 
Chosh (1980); however, our main results differ from theirs. The second 
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related paper Is by Karlin and Rlnott (1980). They introduce a negative 
dependence notion which is closely related to one of ours and they obtain 
some results which are similar to ours. Some remarks about the relationship 
between these two works and the present paper will be given throughout the 
text. 

The main motivation for our definitions is to try to formulate the 
intuitive requirement that if a set of negatively dependent random 
variables is split into two subsets in some manner then one Rubset will 
tend to be 'large' when the other subset is 'small' and vice versa. In 
Section 2 we define the conditions to be discusser}. We derive some in¬ 
equalities in Section 3 and prove some preservation properties in Section 4. 
Examples are given in Section 5. 

In the following "increasing" stands for "nondecreasing" and "decreasing" 
for "nonincreasing". Vectors in are denoted by t * (t,,...,t ) and 

.t <. t_' means ^ <_ t|, i * l,...,n. Similarly t^ < _t' means tj < tj, 
i ■ l,...,n, and 0 *> (0,...,0). A real function on ® n will be called 
increasing if it is increasing in each variable when the other variables 
are held fixed. 

A rv X is said to be stochastically smaller than the rv Y (denoted 

by X Y) if P(X > x) < P(Y > x) for every real x. The random vector 

X ■ (X..X ) is said to be stochastically smaller than Y ■ (Y,,...,Y ) 

— in — in 

(denoted by 3( V Y] if g(X) V g(Y) for every g e C where C is the 
class of Borel measurable increasing functions on 3R n . It' X and Y have 
the same df then we write J( » C Y. It is well known that X \ Y if and 
only if 

P(XeU) £ P(YtU) for every upper Borel set U in Tt n . 


( 1 . 1 ) 









[U is an upper set if x e U and x 5. i implies that x e U.] According 
to Kamae, Krengel and O'Brien (1977), we need only consider open upper 
sets U in JR n . It is also well known that for every random vector X^ 

(1.2) X + a ^ X whenever la ^ 
and that 

(1.3) X + A ^ X whenever A is a nonnegative random vector. 

Also if P(X 0) - 1 then 

(1.4) aX V X whenever a >_ 1. 

See Arias and Lehtonen (1978) for an excellent review on stochastic 
ordering. 

If X and Y are random vectors such that X Riven that Y * y_ 
is stochastically smaller (larger) than X given that Y ■ jp', whenever 
X .1 x' and X and i' are in the support of Y, then we write 

{X|Y - y] 8 t C (?$ X* 

More precisely, this means that for every upper set U, there exists a 
version of P(X^ e u|Y) - <j> (Y) such that $(x) is increasing (decreasing) 
in x on the support of Y. 

A bivariate function K(*,*) which is defined on x (where 
and S 2 are subsets of H) is said to be totally positive of order 2 
(TP^) on x s 2 if K(x,y) >_ 0 and if 

K(x,y) K(x',y') >_K(x,y') K(x’,y) whenever x <_ x 1 , y <_ y' 


(1.5) 
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(see Karlin (1968)). The function K Is said to be reverse regular 
of order 2 (RR 2 ) on x s,, if K(x,y) >_ 0 and If 

(1.6) K(x,y) K(x',y') <_K(x,y') K(x',y) whenever x£x', y <_y', 

(see Karlin (1968), p. 12). 


2. Negative dependence concepts . 

Most, but not all of the positive dependence concepts discussed in the 
literature have negative dependence analogs that can be obtained by changing 
the direction of the monotonicity or of the inequalities which define them. 
Here we define four conditions of this type which we find useful because we 
have methods of identifying distributions which satisfy them. Two of them 
are direct analogs of "conditionally increasing in sequence" and "positive 
orthant dependent". The third has a positive dependent counterpart, but it 
does not seem to have been discussed in the literature, while the fourth is 
a variation of "totally positive of order two (TPj) in pairs". 

In the case of positive dependence, one of the strongest and most use¬ 
ful notions is that of TPj-ness in pairs; that is, the joint density or the 
discrete probability function f is assumed to exist and be TP^ in pairs 
(see Barlow and Proschan (1975), p. 149). The natural negative dependence 
analog then is to assume that f is RR 2 in pairs. There are several 
drawbacks to this notion, however. Firstly, unlike the situation in which 
the joint density (or discrete probability function) is TP 2 in pairs, the 
marginal densities do not necessarily enjoy the same property. A simple 
3x2*2 discrete example suffices to show this. In fact. Theorem 5.1, 
p. 123 of Karlin (1968) ia false when TP 2 is replaced by RRj. (It 
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should also he remarked here that even In the TP^ case, one must make 
some assumptions on the nature of the set {f > 0} before one can use 
the above Theorem 5.1 to conclude that the marginal densities are also 
TP^ in pairs: see Kemperman (1977)). A possible alternative then is to 
assume that not only f, but all its marginal densities are RR^ in pairs. 
But even under this assumption we have not been able to show that our 
weakest condition (see Definition 2.3) is consequently satisfied. Ebrahlmi 
and Qiosh (1980) claim to have proven this result; however, their proof is 
based on an implication which, as will be shown below, does not hold (see 
discussion after Definition 2.4). 

Because of these drawbacks and since we do not always want to assume 

the existence" of a density, we prefer to work directly with the measure 

itself. This point of view is consistent with our other definitions in 

that we do not make any assumptions on the existence of a density in 

defining them. We are thus led to our first definition. 

Let p be a probability measure on the Borel sets in IR n . If 

I,,..., I are intervals in IR^, we define the set function p(I, ,...,1 ) 

in in 

by u(I,,...,I ) ■ p(I x.. ,xi ). By abuse of notation, we write p instead 
in In 

of p. If I and J are intervals in we write I < J if x e I, 

y c J implies x < y, that is, I lies to the left of J. 

2 

Definition 2.1 . Let p be a probability measure on F . We say that p 
is reverse regular of order two (RR 2 ) if 

(2.1) P(I 1 ,I 2 ) Wttj.lp < p(I 1 ,Ip utt{,I 2 ) 

for all intervals 1^ < Ij, I 2 < I^ in F*. We also say that ^(I^.Ij) 
is RR 2 in the variables Ij,!^. If V is a probability measure on F n 












(n >_ 2), we say that y is RR ^, in pairs if y(Ij,...,I n ) is RR 2 in 
the pairs 1^,1^ f° r 1 <_ i < j <_ n when the remaining variables are 

held fixed. The random variables T,,...,T (or the random vector T or 

its distribution function F) are said to be RR^ . in pairs if its 
corresponding probability measure on TR n is. 


Remarks : 

(i) The obvious TP 2 definitions for y are obtained by reversing 
the inequality in (2.1). 

(ii) Clearly, if y is RR 2 (TP 2 ) in pairs, then so are all marginals. 
Furthermore, it is not difficult to show that if F is the 
distribution function associated with y and if 
F(t lt ...,t n ) - yU^,®), (t 2 ,»),..., tt n ,®)) is the survival 
function, then the functions F and F are RR 2 (TP 2 ) in pairs 
in the sense of (1.5) and (1.6). 

(iii) It is easy to show by a simple limiting argument that if y is 
RR 2 (TP 2 ) in pairs, and if y has a density f with respect to 
a product measure m ■ m^*...*!!^ ot o-finite measures such that 
f is continuous on the support of m and zero off the support 
of m, then f is RR 2 (TP 2 ) in pairs. 

(iv) In the n * 2 case, we have the stronger converse, namely, if 

y has a density f with respect to a product measure m ■ m^ x m 
of o-finite measures which is RR 2 (TP 2 ) on * S 2 , where 
is the support of m^(i ■ 1,2), then y is RR 2 (TP 2 ), 

(v) In the TP 2 case, one can generalize to higher dimensions if one 
makes some assumption on the set {f > 0}. Let y have a 
density f with respect to a product measure m ■ mj*...xm n of 
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o-finite measures. Let be the support of m^. Then 

the support of m is S = S . We assume that there 

In 

exists S = S.x...xS such that (f > 0}f\S = S and that f 
In ' 

is TP^ in pairs on S. Then p is TP 2 in pairs. Just 

use Theorem 5.1, page 123 of Karlin (1968) repeatedly to show, 

e.g., that for fixed intervals in 1 r\ 

J n 

g(x lf x 2 ) - /.../ f(x 1 ,x 2 ,x 3 ,...,x n ) dm 3 (x 3 )...dm n (x n > 

I, I 
3 n 

is TP 2 in x^ and x 2 on x § 2 . The result then follows 
by a simple integration using the TP 2 ‘ inequality for g. 

(vi) The generalization to the RR 2 case is not as simple. If one 
assumes, however, that p has a density f with respect to a 
product measure m *> m^x.^xm^ of o-finite measures such that 
the density f when integrated over any n-2 intervals in JR 1 
is RR 2 in the remaining unintegrated variables, then p is 
RR 2 in pairs. In terms of random variables, this can be para¬ 
phrased as follows. Let T^.T n be random variables with a 

density f (with respect to a product measure of o-finite 
measures). Then p is RR 2 in pairs if and only if for every 
1 <_ i < j <_ n the conditional density of 


(T.,T )| r\ (T. e I.} 

1 J k*i,j 

is RR 2 in t^ and tj for all choices of intervals I^(k^i,j) 
in 1R*. Equivalently, if X^ denotes the indicator function of 
I, then p is RR 2 in pairs if and only if 


J.../ [ n 

1R 1R Wi,j 


k 


...t )[ n dt.] 

n k*i,j 
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is RR^ in the unintegrated variables t^ and t^ for all 
choices of intervals I, (k^i,j) in 1R^. By replacing X T 
by <(>^ in the above integral, and requiring it to be RR 2 in 
t^ and tj whenever is a set of functions, 

one obtains the negative dependence condition of Karlin and 
Rinott (1980). It is a condition which is stronger than the 
one of Definition 2.1 as can be easily seen by recalling that 
the indicator function of an interval is PFj. 

(vii) Clearly, if {y^} is a sequence of RRj in pairs probability 
measures and if p n converges weakly to p then p is RR^ 
in pairs. 


Definition 2.2 . The rv's T^,...,T n (or the random vector T or its 

df) are said to be conditionally decreasing in sequence (CDS) if, for 
i = 1,2,...,n-l, 

(2.2) [Ti + i|Ti * t^,...,Tj ■ t^J ♦ (t^,...,t^). 

Definition 2.3. The rv's T,,...,T (or the random vector T or its 

- i n — 

df) are said to be negatively upper orthant dependent (NUOD) if for every £, 

(2.3) P(T > t) < " P(T, > t,). 

— — — i-l i i 

They are said to be negatively lower orthant dependent (NLOD) if for every t^, 

P(T < t) < it P(T. < t.). 

-- i-l i - i 


(2.3') 
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When n * 2, (2.3) and (2.3') are equivalent, but not when n >_ 3 
(see, for example, Ebrahimi and Ghosh (1980)). 

The next concept has a natural positive dependence analog, however, 
we are not aware of any place in the literature in which it has been 
discussed. 


Definition 2.4 . The rv's T^,...,T n (or the random vector T. or its 
df) are said to be negatively dependent in sequence (NDS) if, for 
i = 2,3,... ,n, 

(2.A) [(Tj,...,T ± _ 1 )|T ± - t ± ] ti . 

Often, to verify (2.4), one can find it easier to verify that for 
i ~ 1 ,...,n, 


(2.4') 


I »• • • »• • • I 


t i 1 


i- 


Then clearly (2.4) holds. 

We now investigate some of the relationships among the various 
definitions. First note that NDS implies both NUOD and NLOD and these 
implications are sharp. To see this, use methods similar to Barlow and 
Proschan (1975) to show that (2.4) implies for 1 ■ l,...,n-l. 


(2.5) 


[(t 1 ,...,t 1 )|t 1+1 > t i+1 


] ? [(t x .Tj) |T 


i+1 * t i+l 1 


whenever tj + ^. (Although Barlow and Proschan assumed the existence 

of a density, a modification of their proof works.) But from (2.5) it 


follows that 
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P (T. > t 


1*' 




C n> 


i P(Tj 


.,T 


n-1 


t P(T n 

n-1 n 




i P(Tj > tj 


••V2 » ‘„-2> 


J P(T » t,) 
i-n-1 


i ... i " P(T. > t.) 
1-1 1 1 


which proves (2.3). The proof of (2.4)=>(2.3’) Is similar. 

It is not difficult to construct an example which shows that 
CDS NUOD and that CDS NLOD and hence CDS^NDS. For example, let 
P(T 1 -1,T 2 -1) - P(T 1 -l,T 2 -2) - P(T 1 -2,T 2 -2) - .1 and PCT^.Tj-l) - .7 
and let given T^ - 1, T 2 - 1 be degenerate at 11, given 

Tj - 1, T 2 -= 2 be degenerate at 1, T 3 given Tj - 2, T 2 “ 1 be 
degenerate at 10 and T 3 given T^ - 2, T 2 - 2 be degenerate at 1. 

Then, clearly, (T^T^T-j) is CDS but P(T > 1, T 3 > 1) - .7 > .64 
= P(T ] > 1) P(T 3 > 1), thus (T 1 ,T 2 ,T 3 ) is neither NUOD nor NLOD. 

Fbrahimi and Ghosh (1980) claim that CDS=>NU0D; the example shows that 
this is not the case. 

Next we show that NDSjt^CDS. Let (T^,T 2 ,T 3 ) take on values on 
the eight vertices of the unit cube such that 

PfTj-l, T 2 -0, T 3 =l) - KTj-1, T 2 -l, T 3 -0) - .2 and the other six pro¬ 
babilities are .1. It is easy to see that and T 2 are independent 

and that [(Tj.T^Tj - 0] V [(T^T^Tj - 1]. Thus, (T^Tj,^) is 
NOS. But P(T 3 > 0|T x - 0, T 2 - 0) - 1/2 < 2/3 - P(T 3 > 0|T 1 - 1, T ? - 0), 
hence (T^.T^T^ is not CDS. 

We will shortly show that under some reasonable assumption, RR 2 :=£CDS. 
It is not known whether RR 2 ^ NDS, but it does imply NUOD and NLOD. This 
is clear since from the remark (ii) following Definition 2.1, we have that 
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hot!) K and F are RR 2 In pairs. The result then follows by a 
simple argument. 

The implications RR 2 =£CDS, RR 2 =?NUOD, RR^ NLOD, NDS=£NUOD 
and NDS^NLOD justify the consideration of the RR 2 and the NDS 
concepts. In Sections 4 and 5 it will be shown that many df's are RR 2 
or NDS and thus these df's satisfy the meaningful CDS concept and the 
inequalities which can be derived from the NUOD and NLOD concepts. In 
addition to it the NDS concept is intuitively meaningful by itself. We 
mention, in passing, that since Karlin and Rinott's (1980) condition 
implies the RR 2 condition it follows that it implies the NUOD and the 
NLOD inequalities (2.3) and (2.3'). Actually, Karlin and Rinott (1980) 
obtained some additional useful inequalities which follow from their 
stronger condition. 

To justify calling (2.1) — (2.5) "conditions for negative dependence" 
we have to show that they imply 

(2.6) COVOl^.Tj) <_ 0, i £ i < J £ n, 

when the second moments exist. 

From (2.3), it follows that P(T^ > t^, T^ > t ^ ) £P(T^ > t^) P(Tj > tj) 
and it is well known that this inequality implies (2.6) [see e.g., Lehmann 
(1966)]. Similarly (2.3*) implies (2.6) and hence also (2.4) and (2.5) 
imply (2,6). 

If 0^,...,^) is CDS then [T^^ - tj hence covO^,^) <0. 

Thus, if (T ff (i),... ,T^ (n)) is CDS for every permutation tt of {l,2,...,n} 
then (2.6) holds. 

Vie close this section with a proof that RR 2 in pairs implies CDS under 
some reasonable assumptions. Let p be RR 2 in pairs and let (T^,...,T^) 





1 'f anv random vector having 11 as its inducer) probability measure. 


lix t l IK and an integer i, 1 <_ i n-1. Define the set functions 
v and X on the Borel subsets of TR* by 

v(A) - P(T 1+1 > t 1+1 ; (T x .T t ) e A) 


A (A) - P((T .Tj) e A). 

£■ °° 

Note that X is just a marginal of y. Then if is a sequence 

of rectangles which partition TR* for each l and whose mesh size tends 

to zero as £-*■<*>, and if (R,, } is a refinement of (R,, } for all £, 

ik ik 

then it is well known by martingale arguments (see, e.g., Meyer (1966)) 
that on a set I) with X(D C ) * 0, v ^ik^^ R ik^ "* 38 

decreases to (t^.t^) pointwise and in L^; moreover <J>(tj,...,t^) 

is a version of P < T 1+1 > t i+il T i * c i. T i * t^)* We may assume without 

loss of generality that D C support X. We want to show that <p is 
decreasing on D. This will follow if we can assume that the support 
A of X satisfies a chain condition; that is, if t e A and t/ e A witli 


t < t', then there exist t, e A such that t = t^ < t. <...< t “ t' and 

— — — —j — —i) —1 — — —m — 

t, differs from t... in only one component. To see this, suppose that 
“1 “ 1 + 1 

t^ and _t differ in only the first component. Then let 

h ‘ h+1 ■ ‘-•'W * S + 1 ■ (t i+i'"> and h . h be "v 

intervals such that e 1^ * —1 E ^l x l2 X, '* x *i‘ p 

is the marginal of (Tj.T j), then since it is RR 2 in pairs, we have 











.• '■■ 
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or that 

v(Ij *^ 2 * * * * ^^1 *^"2 *''' '^1 ^ 

uij.ij. 1 xa'.ij,...,^) • 

Note that there is no trouble dividing since the denominators are non zero. 

Iterating this procedure to pass from t. to t_, t- to t, ,...,t . to 

t and then letting the intervals shrink we find that <fr(t) > 4>(t'). Now 
~m — — — 

set 

C$(t) if t^ e D 

<(,*( 0=1 

inf{$(s) ; £ e D, s_ <_ _t} (inf 0 - +=») if t^ D. 


It easily follows that 4>* is decreasing everywhere and is a version of 
P * T i+l > C i+1I T 1 “ t l’*“’ T i = t i^* 


The chain condition is easily seen to be satisfied if the support 
of n is a cross product. Without some type of chain condition, we can 
only show that <(> is decreasing componentwise on D. In this case, one 
may not be able to extend <(, such that it is decreasing everywhere and is 
still a version of the condition probability. 


3. Some inequalities . 

This section is devoted to the derivation of some inequalities which 
may be of special interest. The results closely parallel those found in 
Karlin and Rinott (1980), but are derived under weaker assumptions. 

Proposition 3.1 . 0^.T^) is NUOD if and only if 

E{ " *.(T.)] 1 " E[* (T )] 

i=l 11 i-1 




(3.1) 
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wlienever all $ are nonnegative and Increasing. The result (3.1) 
is also true if we replace NUOD by NLOD and increasing by decreasing. 

Proof. Let (t) - X /t v (t), i - l,...,n. Then (3.1) reduces to 

- i (bjt”) 

the NUOD inequality (2.3). Since each side of (3.1) is multilinear in 
the the result holds for nonnegative linear combinations of such 

indicator functions and hence for the general by a standard limiting 

argument. | | 

Now suppose that p is RR 2 in pairs. For every i, let 1^ ■ 
all intervals, with < K^. 

Theorem 3.1 . If 1 5_ k <, n, then 

0.2) uOp.-.v .(li.—.v i.Uj.\.i k+1 .V 

X U<I 1.*k' J k+1.V - 

The result (3.2) is also true if ve replace all J's by K's. 

Proof . We proceed by induction. If n * 2, then 

fliU.,!,) ii(I,,JJ if k-1. 
u(J lt J 2 ) wd^ij) i) 1 

|_p(J lt J 2 ) p(I 1 ,I 2 ) if k - 2. 

The case k ■ 1 follows from the RR 2 assumption since 

u(J 1( J 2 ) (Ji.J 2 V.K 2 ) 

» 

(J^^.Jj) (JjyKj.J^Kj) 


0 > 


|(J 1 ,J 2 ) y(J lt K 2 ) 


VJ (K x ,J 2 ) p^.Kj) 










and the case k = 2 is an identity. 


Now suppose that (3.2) is true whenever v is a probability measure 

on lR n which is RR 2 in pairs. Let p be a probability measure bn lR n+1 

which is RR^ i*i pairs and let 1 <_ k <_ n+1. Since there is nothing to 

prove if k = n + 1, we may assume that 1 <_ k <_ n. Similary, we may 

assume that p(J,,...,i ,J .,) 4 0. It suffices then to prove that 
l n n+i 


(3>3) - (I 1.VW _ - (I 1.vv . 1 _ 

" (J 1. J k ,I lctl.V'rrt-l* “ P<J 1. J k’ I k+l. I n’ J irt-l ) 

since by the induction hypothesis, we have that 

ti(I l. I n ,J n+l ) _ < tl(I l. T k ,J k+l. J n ,J n-H ) 

p(J 1 ,...,J k ,i k+1 ,...,r n ,J n+1 ) - w(J 1 ,...,J n ,J n+1 ) 

But, 

u(I l , ” ,,I n ,I n4-l ) y(J l ,I 2. I n ,I n-H ) 

u(I 1 ,...,I n ,J n+1 ) ~ udidj. I n ,J n+l ) 

p(J l. J k ,I k+l ,,,,,I n ,I n+l ) 

p(J 1 ,...,J k ,I k+1 ,...,I n ,J n+1 ) 

which is another way of writing (3.3). The jth inequality above follows 
from the fact that p is RR 2 in the pair j and (n+1). || 

Suppose that (T ,...,T n ) is RR 2 in pairs. 

Corollary 1 . If a, 6 partition (l,...,n), then 
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(3.4) P(T i eJ 1 ,iea^e) P(T^el t ,ica’u B) 

< P(T i eJ i ,iea;T j eI j ,j€6) PO^e^,ie«;T cJ^,jeB). 

It also holds true If we replace all J's with the K's. 


Corollary 2 . If a, 6, y partition {l,...,n}, then 

(3.5) P(T 1 eL i ,iea;T j eJ j ,jEBuY) P(T 1 eL i> iea;T j eI j .jeBoy) 

^ P(T^eL^ ,lea;Tje.J eBjT^eI^, kEY) P(T^eL^, iea;TjElj, jeB;T k EJ k ,kEY) 

for any intervals L^, i e a. It also holds true if we replace all J's 
with the K's. 

If we take ° (-“.bjJ and 1^ - (-“>,“>) for i £ au B in 

Corollary 1, we get (3.5) of Karlin and Rinott (1980). If we take, in 
Corollary 2, = [a^.b^], i e a, ■ (-“.b^] and 1^ ■ (-®,®) for 

1 e Rvj y w e get (1.7) of Karlin and Rinott (1980). 

Also, note that as soon as we have an inequality of the form 


<3 - 6) P(T 1 i b l. T „ i b „> i P(T 1 i b l. b k> P(T k« iVl. T „ i b n> 

it follows as in the proof of Proposition 3.1 that 


(3.7) 


E[ tt <*.. (T.)] < E[ ir <MT.)] E[ tt d> (T.) ] 
i-1 1 1 ” i-1 1 1 j-k+1 i J 


whenever 4^ are nonnegative and decreasing. Similarly, if we have 


P(T 1 » a l. T „ » *„> 2- r < T l * »1. T k * V P(T k« * “k+1. T „ * *n 


(3.8) 





then (3.7) holds for $ nonnegative and increasing. 

In particular, (3.6) holds If F is RR,, In pairs and (3.8) holds 
if F is RR£ in pairs. In fact, if F(or F) is RR,, in pairs, then 
it Is MRR^ in the sense of Karlin and Rinott (1980); i.e., 

F (* A Z.) F(xvy_j0 <F(x) F( Z ). 


4. Closure Results . 

Preservation theorems are useful for identifying negatively dependent 
df’s or for constructing new negatively dependent df's from known ones. 

In this section we discuss some preservation results and describe a method 
for the construction of negatively dependent df's. 

Theorem 4.1. If T,,... ,T are (*) and if 1 li, ,... ,\p are strictly 
in in 

increasing functions then ^-(T.),...,♦ (T ) are (*) where (*) is 

lx n n 

one of the following: RR 2 in pairs, CDS, NDS, NUOD or NLOD. 

Theorem 4.2 . If (T^,...,T n ) and (S^,...,S n ) are independent and are (*) 
then (T^ T^, S^,...,S n ) is (*) where (*) is the same as in Theorem 4.1 

The proofs of these theorems are straightforward and will be omitted. 
The following preliminaries are needed for the statement of Theorem 4.3 

A univariate density f is said to be a Polya frequency function of 
order 2 (PF 2 ) if f(x-y) is TP 2 on H * ®. A probability function f 

is PF 2 if f(x-y) is TP 2 on f! x IN where V m {...,-1,0,1,...}. A 


thorough discussion of PFj densities and many examples can be found in 
Karlin (1968). 
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It will ho shown In a later section that many multivariate random 
vari.il les satisfy a certain structural condition. In the next theorem 
we will show that this condition implies some of the dependence conditions 
that we have introduced. 


Theorem 4.3 . Let Sg,S^,...,S n be independent rv's and assume that 

each has a PF 2 density (or probability function). Fix s and let 

(T, ,.. . ,T ) have the same joint df as the conditional df of (S,,...,S ) 
L n in 

given that S. + S,+...+S « s, that is, 

() 1 n 


(4.1) 


(T 1 ,...,T n ) S - I(S 1 ,...,S n )|S 0 + S 1 +...+S n - s]. 


Then (Tj,...,T ) is RR 2 in pairs and consequently CDS, NUOD and 


NLOD. 


Proof, Let y be the probability measure of (T^.T ) on lR n . Then 

by assumption p(I,,...,I ) • P(S, e I,,...,S e I IS_ + S-+...+S ■ s). 

L n ii nnui n 

Now the loint density of [(S,,...,S )IS rt + S-+...+S ■ s] is given by 

i n u i n 

n 

C TT f (s )f (s - S--...-S ) 
i= l i i 0 1 n 

wliere c is a normalizing constant. We first show that p is RR 2 in 
the variables Ij, I 2 when the remaining intervals I^,...,I n are held 
fixed. According to the remark (vi) following Definition 2.1, we need 
only show that 

g(s 1 ,s 2 ) - c f^s^ f 2 ( s 2 ) f o^ s_s l’ ,8 2^ 


is RR 2 in s^ and s 2 » where 
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• '('.) = /•••/ f.(s )...f (s ) f (f-s _-...-s ) dm(s )...dm(s ) 
<> 1 3 3 nn03 n 3 n 


r T T 

) n 


and m is either the Lebesgue measure or the counting measure. However, 

t lie above is nothing but the convolution of the PF^ functions 

f, X. .f X. and f_, where X is the indicator function of the 

i I n I u A 

3 n 

set A, and so f^ is eas *-ly follows then that g is R R 2* 

The proof that u is RI^ in the variables 1^, 1^ for all 1 £ i < j n 
is similar. I| 


Remark . Actually, it is not difficult to show that under the assumptions 

of Theorem 4.3 the df of the random vector (T ,... ,T ) satisfies the 

1 n 

stronger condition of Karlin and Rinott (1980). The consequences of this 
ole ervation are discussed below. 

It is well known (see Section 5) that the multinomial, the multivariate 
hypergeometric, the Oirichlet and the Dirichlet compound multinomial random 
vectors as well as some multivariate normal random vectors with nonpositive 
correlations, can be represented as in (4.1). Thus, all the respective df's 
are RR^ in pairs. Karlin and Rinott (1980) have shown that these df's 
actually satisfy their stronger condition; however, their proofs are quite 
involved and differ from one case to another. By the Remark after Theorem 4.3 
these results of Karlin and Rinott (1980) follow at once. 

The following theorem will be found useful in the next section. 


Theorem 4.4 . Assume that (T^,...,T n ) and (S^,...,S n ) are independent 
and NOS. If all the univariate marginal densities (with respect to Lebesgue 
measure) or probability functions in the discrete case of S_ and T are 
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PF„, then (T, + S. ,...,T + S ) is NDS. 

2 1 I n n 


Remark . Karlin and Rinott (1980) have proven a similar result. They 


assumed that IS and T satisfy their RR 2 condition and that they 
have PF^ marginals, and they showed that then £ + T satisfy some 


inequalities that are essentially variants of the NUOD and the NLOD 


inequalities. 


For the proof of Theorem 4.4 we need the following lemmas. 


Lemma 4.1. Let X and Y be independent rv's and assume that Y has 


a PF^ density (probability function). Then 


[X|X + Y - z] z. 


This lemma is actually Example 12 of Lehmann (1966). 


Lemma 4.2. Let X and Y be independent and assume that X and Y 


have PF 2 densities (probability functions). Then 


(4.2) 


[(X,Y)|X + Y - z] z. 


Proof. Let z-^ < z n . Denote by U (x,y) the upper set {(s,t): s > x, t > y). 


First it will be shown that for every (x,y) 


,3) P{ (X ,Y) e U (x y) |x + Y * z i } L p UX,Y) e U (x y) |x + Y - z 2 ). 


Let “ {(s,t): s + t - , x < s < - y}, 

B - {(s,t): s + t - s > z ± - y), C ± - {(s,t): s + t - z ±t s < x). 


I 


I 







i = 1,2 


I s c 

(see Figure 1). From Lemma 4.1 [X|X + Y = z] t z; hence 


P{ (X,Y) e C 1 |X + Y = z 2 ) = P{X <_ x|x + Y = z } 


>. P{X < x|x + Y - z 2 ) - P{ (X,Y) e C^jX + Y ” z 2 >. 

Similarly 

P{ (X ,Y) e B^X + Y - z x } > P{(X,Y) t B,,|x + Y « z 2 >. 

Tlius 

P{ (X,Y) e A X |X + Y - Zj) > P((X,Y) e A 2 jx + Y - 
which is (4.3). 

t 



Figure 1 . 

m 

Now consider upper sets of the form U * W U, \ * called funda- 

i«l 

mental upper domains in Block and Savits (1979). Without loss of generality 
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assume that x^ < <_ .. • x m and Yj i i 


> y . Define 
— m 


•* min {i: x^ + y <_ z^} 

* min {i > i^: x 1 + y i > z^s 

and, by induction 

(A.4) i R+1 - min {i > j R : x t + y ± <_ z 2 > 


( 4 . 5 ) 


k+1 


min {i > i 


k+1' 


X i + y i > Z 2 } 


k = l,2,...,n where n < °° is the largest k such that the set on the 
right hand side (RHS) of (4.4) or (4.5) is not empty. If there are n i's 
but only n-1 j's define “ m + 1 (see Figure 2 in which m ■ 12, 

n - 3, i x - 3, Jj - 6. i 2 " 8, j 2 - 11 , i 3 - 12 , j 3 - 13). 

n 


from 


Let U « U ..In the following the first inequality follows 

lc»l ’ 

*k 3 k-l 

UC U and the second one from (4.3): 


P{(X,Y) e U|x + Y - z L } < P{ (X,Y) e U|X + Y - z^ 


n 

- I P{ (X,Y) e u Jx + Y - z } 

k ‘i ( W 1 

(4.6) 

< l P{(X,Y) e U JX + Y - z,} 

k-i W 1 


P{ (X,Y) e U|X + Y - z 2 ) - P{(X,Y) e u|x + Y - z 2 ). 












Si net* open upper sets can be approximated by fundamental upper domains (Block 
and S.ivits (1979)) it follows from (4.6) that for every open upper set 1' 


P{ (X,Y) c U|X + Y - Zl > <_ P{ (X,Y) c U|X + Y - z?) 
and the proof of the lemma is complete by (1.1). || 

Lemma 4.3. Let X * (X, ,...,X_) and Y = (Y,,...,Y ) be independent 
--— — I m — l m 

and assume 


(4.7) 


f(XX i)IX - x ] V x 
x m-x m in i 


and 


(4.8) 


[(Y ...,Y ,)|Y - y ] V y 

i m—l mm m 


Furthermore, assume that X and Y have PF_ densities (probability 

mm / 

functions). Then 


[ (X + Y .X + Y ) | X + Y - z ] z, 

i x m—x m—± m m m i 


Proof . Dearly, for any increasing function g, 

F(g(X + Y ? X + Y ) |X + Y - z ] - F[<KX ,Y )|X + Y « z ] 

1 I m—x m-x m m m mm'mmin 

v/here rp (x ,y ) «= Ff g(X, +Y,,...,X . + Y .)|x **x,Y - y ]. However, 

mm 1 1 m-x m—x m mm m 

■P(x ,y ) decreases in x and in y because of (4.7), (4.8) and in- 
m m m m 

dependence. Thus, by Lemma 4.2, E [ (X ,Y )|X + Y » z ] decreases in 

m m • m m m 


7 . . 

m 


Proof of Theorem 4.4. Let i e {l,...,n-l}. Substitute m ■ i + 1 in 
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! omnia 4. 1 to obtain 


[O', + S, 


• ,T. + 


S i ) 


i+1 


+ s 


i+] 


Z i + 1 ! 




i+l ’ 


that is, T + S Is NDS. |j 


5. Examples . 


r >. 1 . 'Hie multinomial df . 

l-et (Tj,...,T ) have the loint prohahility function with parameters 

(^iPj i•••») j 


P(T 1 ■ V 


,T = t ) = 


n: 


n n 


t 1 :,,,t n : ( N ’j 1 t 1 ) : 


n 

n- ( iii p « > 


i-1 
n 


N " l 

n 1 n 

* ( 1_ l Pi> . t > 0, [ t < N, 

1-1 1-1 1 


where p > 0 (1 • l,...,n) and 0 < £ p < 1. 

1 1-1 1 

The multinomial df Is the conditional df of independent Poisson 

rv’s Riven their sum. Thus, by Theorem 4,3 the multinomial df Is RR^ 

in pairs and hence it is also CDS, NUOD and NLOD. By Remark (ill) the 

joint probability function of (T ,... ,T ) is RR„ in pairs. Ry the 

in L 

discussion after Theorem 4.3 the multinomial df satisfies the RR^ 

condition of Karlin and Rinott (1980). 

To show that (T,,...,T ) is NDS it is enough to show that for 
1 n 

n < t < N-l, 

- n — 


( 5 . 1 . 1 ) 


f (T, 


T ) 
’ n-l' 




f (T, 


r .)|t 

n-l 1 n 


t+1 ] 
n 








mse then, l>y symmetry one has (2.4'). The left hand side (l.HS) of 

(5.1.1) has an (n-1)-dimensional multinomial df with parameters 

(N-i n ,»| i ,. . . ,q j ) and the right hand side (RHS) of (5.1.1) has the 

same df with parameters (N-t -1, q,,...,q .) where q, ■ p./(l-p ). 

n i n-i 11 n 


Thus, 


(5.1.2) 


st 


LHS (5.1.1) - RHS (5.1.1) + (S,,...,S ,) 

i n-l 


where (S ,...,S .) has a multinomial df with parameters (1,q ,...q 

l n-l 1 n 

and the sum on the RHS of (5.1.2) Is of independent random vectors. 
Inequality (5.1.1) follows now from (5.1.2) and (1.3). 


l.et (X^,...,X ) he th? sum of Independent n-dimensional, not 
necessarily Identically distributed, multinomial random vectors, that is 
let (X^,...,X^) have the probability generating function 


m N * 

<P u u i + — + P n * U n ) 


n 


where p^ ^ 0, (i = 1 

,...,n) 

, 0 < 1 P U < 1 
i-l 

and >_ 1 (4*1,2,. 

The univariate marginal 

df's 

of a multinomial 

df is a binomial df. 

Since the binomial df 

has a 

PF 2 probability 

function it follows from 

Theorem 4.4 that (X^ ,. 

-V 

is NDS. Hence 

(Xj,...,X n ) is also NUOD 


and NI.OI), that is, 


and 


P(X L > x lf ...,X n > x n ) < P(X i > x t ) 

p < x i — x i.i*„> i,S, p «,iV- 


These inequalities are stated unproven in Mallows (1968). Compare also 
l.ehmnnn (1966), pp. 1 143, 1144 and 1151. 
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•. ■'. Mu 1t i v.-»r ia t e normal. 

lot T = (T,.T ) lip ;; multivariate svnmerrie normal ramlom 

1 n 

vet tor with C.orr (1\ ,T ) * p < 0, 1 <_ < i * n. Then r> > - (n-1) 

Ive will show that T is RR 2 in pairs. 

Using Theorem 4.1 assume, without loss of generality, that FT^ = 0 

anti Var(T ) = 1, 1 = l,...,n. I.et Y,,...,Y he independent identically 
1 in 

distributed normal rv's such that EY^ = 0 and Var(Y^) « 1-p (i - 1.n) 

and let he an independent normal rv with F.Y^ - 0 and 

Var(Y 0 ) = (-p)" T (l-p)(l+(n-l)p). Then 


< T i--V - [(Y 1 


,Y n )|Y 0+ Y] + ... + Y n 


0]. 


Since any normal density is PF,, it follows, from Theorem 4.3, that T 
is UR,, In pairs. 

In fact we can obtain a stronger result. Tf the correlation matrix 
of T is of the form 


(5.:\i) 



n \ 


^(r r l) 1/2 \ 

V 

r ni 


Uv 1 > 1/2 I 


,, , xl/2 ,,1/2, 

((r^ ~3 ) «••• » 1) ) 


wIkti* 


pi i rs. 


I ,... , n a nd 


n 


l 

i-1 



> n-1 


then T Is RR^ in 


To show it note that every matrix of the form (5,2.1) can he the 

correlation matrix of a multinomial random vector, X_ = (X^,,,.,X ), sav. 

(?) 

bet X_ , {, = 1,2,..., be a sequence of independent random vectors 

( N ^ / Q \ 

distributed as X_. Clearly Y/ ' * J It is a multinomial random 

i- 1 , * 

vector with correlation matrix (5.2.1). Normalizing Y' such that it 






has zero moans and unit variances, it converges in distribution, hv the 
hi I t i v.ir i ate central limit theorem, to a multivariate norma] random 
vector with correlation matrix (5.2.1). By Remark (vii) the limit in 
distribution of UR^ in pairs random vectors is RR^ in pairs. The 
assertion in the preceeding paragraph now follows. The previous result 
that deals with the symmetric multivariate normal df with negative 
correlations p is obtained by taking r * 1 - p in (5.2.1). 

Every multivariate normal random vector T with nonpositive 
correlations is NDS. It is easy to verify (2.4') directly using (1.2). 
Thus T satisfies (2.3) and (2.3'). 


5.3. M ultivariate hypergeometric . 

let (Tj,...,T ) have the probability function 

/m- y m,\ 

■i 
: i 


l’(T = t 1 ,...,T n - t n ) - VNi 


M. 


n 
it 

i*=l t 


1 hi 

\ i-1 7 


1 i-1 


with positive integer-valued parameter vector (N, M ,... ,M , M) (see 

1 n 

Johnson and Kotz (1969]. 

The multivariate hypergeometric df is the conditional df of 

independent binomial rv's given their sum. Thus, by Theorem 4.3 the 

hypergeometric df is 1^ in pairs and hence it is also EDS, NUOD 

and NLOD. By Remark (Hi) the loint probability function of (T, ,...,T ) 

l n 

is KR^ in pairs. A special case of this fact was observed by Lehmann 
(1966), p. 1144. See also Ebrahimi and Ghosh (1980). By the discussion 
after Theorem 4.3 it follows that this df satisfies the RR^ condition 


of Karlin and Rtnott (1980). 
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We will show now that T fs NI)S hy proving that, for 0 < t < N-l, 

— — n — 


<5.3. 1) 


<T i. T „-i>l'„ ‘ i i< T ,.Vi>! T « ’ * 11 


Then, by symmetry, (2.A 1 ) follows. Since the random vectors on each side 
of (5.3.1) have multivariate hypergeometric df we actually have to show 


(5.3.2) 


U V V 


where U - (U. ,...,U ^) has hypergeometric df with parameters 

(N + l, M ,...,M ,, M + 1) and V ■ (V. ,. . . ,V .) has the same df with 
1 n-I — 1 n-l 

parameters (N, ,, M) where N ■ N - t - 1 and M - M - t - 1. 

1 n-l n n 

Denote m ■ n - 1. 

Thinking about as the number of individuals in the sample in 

Category i (i = 1,...,m) it is easy to see, by conditioning on the category 
of the first individual chosen, that 

(Uj.U m ) S . C (w{ l) ,... t w£{, wj £) + 1, wj*}.W< 4 >) if the first 

s t 

individual is in category £, £ - 1,2,...,m, and that (U,,...,U ) » (V,,...,V 

l m l m 

If tilt* first individual is in neither of categories l,2,..,,m, where 
(t) (l) 

(W ...,W ) has a multivariate hypergeometric df with parameters 

1 m 

(N, j, M-l, M). Thus for proving (5.3.2) it 

suffices to show that 


(5.3.3) (WW<°+1, W. W U)) s^t . V J, £ - 1 .m. 


We will prove (5.3.3) when £ ■ 1; the proof for the other t’s is 
similar. Omit the superscript 1 and consider (W^,...,W m ) ■ (W^\...,W^ ^ 











which has ,i hypergeometrio distribution with parameters 



Then, it is easily seen, that the rv’s »•••» v m can have che following 

interpretation: first divide a population of size M into a group of size 

N and another group of size M - N. Next, choose at random individuals 

and let be the number of them in the first group (of size N), then 

choose at random individuals of the remaining M - individuals and 

let V 2 be the number of them in the first group. Continuing this way, 

finally choose M individuals of the remaining M - M. - ... - M . 

individuals and let V be the number of them in the first group. 

m 

The rv's W,,...,W may have a similar interpretation. 

1 m 

With these interpretations it is easily seen that if the first 

individual chosen in the 'V experiment.' is not in the first group then 

V = l W S < C (W, + 1, W„,... ,W ), where the inequality follows from (1.2). 

— - 1 l m 

otherwise (V t .V m ) (U } + 1, W 2 ,...,W m >. Thus, unconditionally 

v : '< l (Wj + l,...,W m ), which proves (4.3.3) when 9 . - 1. 







5.4. 


The Dirichlet distribution. 


l.et T = (T.,...,T ) have the density 
— In 


r( l ej 


f (tt ) — 

1 n n 


L i' 
1-0 3 


n 8-1 n 9,-1 

o - It ,) 0 . t , J 


* r(o ) 

i-o 3 


l t,) » t 3 , t > 0, l t < 1. 

j-1 3 1=1 1 3 j-1 3 


whore the parameter vector (0^,0^»* • • *6 n ) satisfies 0^ > 1, j = 0,1.n. 

The Dirichlet <lf is the conditional df of independent gamma rv’s 
given their sum. Thus, hy Theorem 4.3 the Dirichlet df is RR 2 in pairs 
and hence it is also CDS, NUOD and NLOD. By Remark (iii) f is RR 2 in 
pairs. Special case of this fact is Example 10 (iii) of Lehmann (1966). See 
also Ehrahimi and Ohosh (1980). By the discussion after Theorem 4.3 it 
follows that this df satisfies the RR 2 condition of Karlin and Rinott 


(1980). 


To prove that T is NDS it is enough to show that for t < t' 
- n — n 


0.4.1) (<T,. T n -l>l T n " ¥ t (T,, 


•Vi>i T n - 


then (2.4') follows. 


It is easily seen that 


(i-t n )-‘i<T 1 ,.... Vl )h n - t„) s -' d-t^)- 1 ((T 1 .- t;i 


Since 1 - t > 1 — t * (5.4.1) follows from (1.4). 

n - n 


5.5. Dirichlet compound multinomial df. 


Let T - (T,,...,T ) have the probability function 
- I n 







..V- X r.’ 


P(T 


i l r 


,t = t ) = 

n n 


N) r '.LV . r(t +b , 


n 

r(N + l o ) 

i=o 1 




1=1 'r r( V 


3? 


r(N - j'j + V 

-f- 

(n - l t.):r(e.) J 


j-i 


where N is a positive integer and 0^ >_ 1, j ■ 0,1,...,n [see Johnson 
and Kotz (1069)]. 

The Dirichlet compound multinomial df is the conditional df of 
independent Pascal (negative binomial) rv’s given their sum. Thus, hv 
Theorem 4.3 this df is RR^ in pairs and hence it is also CDS, NUOD 
and NLOD, and by the remark after Theorem 4.3 it satisfies the RR^ 
condition of Karlin and Rinott (1980). 

To show that T is NOS it suffices to prove that 


(5.5.1) 


[(T ...,T )|T - t ) t 

1 n-i n n n 


then (2.4') follows by symmetry. But the random vector in (5.5.1) has 

the Dirichlet compound multinomial df with parameters N - t^, 

0„,...,() Using this fact it is easy to verify (5.5.1). 

0 n-i 


5.6. Df's supported by negatively tilted surfaces . 

In this subsection we illustrate two simple examples of such df's in 
The general idea then should be clear and will not be detailed here. 
I.et (T] ,^2 ,T^) have the uniform df on the surface 
{(t r t 2 ,t 3 ) 1 (0,0,0): t x + t 2 + t 3 - 1). Then f (Tj ,T 2 )|T 3 - tj] has 








||. 1 r ... ■ - iTffl 


ww^rr^. . pgp m i 
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a uniform df on the segment — (0*0): 1 1 1 2 = 3 - t^) 

which is clearly stochastically decreasing in t^. By symmetry (2.4') 
holds, hence ( T i ,T 2 ,T 3^ is NDS. Clearly ( T i»^2 ,T 3^ is also CDS. 

Similarly, if (T^.T^.T^) have the uniform df on the surface of 
the unit sphere intersected with the positive (or the negative) orthant 
then (Tj.T^.T^) is NDS and CDS. 

5.7. N UOD df's in reliability theory . 

Buchanan and Singpurwalla (1977) considered the class of nonnegative 
multivariate new better than used (NBU) df's that satisfy, for all 

s > 0, £ > 0, 

(5.7.1) PIT, > s , ♦ t,.T, > s n + t n ) < PIT, . T, . .„) 

“ r(T l > h. T „ * '„>• 

It is well known that, when n » 2, (5.7.1) implies that ("’'i*'^) are 

negatively quadrant dependent, that is 


(5.7.2) 


P(T 1 > t l' T 2 * t 2 ) i P(T l > t l ) P(T 2 > t 2 )> 


h'(' will show, by induction, that in general (5.7.1) implies that 

(T,,...,T ) is NUOD. From (5.7.2) we know it is true when n - 2. Assume 
1 n 

that 


P(Tj 


,T 


n-1 


> 


t 

n- 


1 


) 


< 


n-1 

ir P(T > t.). 
i-1 1 1 


(5.7.3) 






Substitute 


Lo obtain 
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